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We show that the S=l/2 pyrochlore lattice with both Heisenberg and antisymmetric, 
Dzyaloshinsky-Moriya (DM) interactions, can order antiferromagnetically into a state with chiral 
symmetry, dictated by the distribution of the DM interactions. The chiral antiferromagnetic state 
is characterized by a small staggered magnetic moment induced by the DM interaction. An external 
magnetic field can also lead to characteristic field-induced ordering patterns, strongly dependent on 
the field direction, and generally separated by a quantum phase transition from the chiral ordered 
phase. The phase diagram at finite temperature is also discussed. 



I. INTRODUCTION 

The behavior of many-body systems involving quan- 
tum spins has been one of the central topics in recent 
years since the properties of such systems are relevant to 
a great variety of materials, mostly oxides. The structure 
of the ground state and the various symmetry broken 
phases that emerge are issues of special interest, espe- 
cially in systems of low-dimensionality and/or where frus- 
tration is present i In this context the Heisenberg model 
on the three-dimensional pyrochlore lattice consisting of 
corner sharing tetrahedra, shown in Fig. ^a), is in a 
league of its own. The pyrochlore lattice is strongly ge- 
ometrically frustrated and is relevant to numerous com- 
pounds. It has been argued that no magnetic order is 
present in the ground statei^i 3 ^ The effects of various ad- 
ditional interactions have also been studied, such as mag- 
netoelastic couplings, 5 long-range dipolar interactions, 6 
and orbital degeneracy^ These interactions (in addition 
to various anisotropies) can generally lead to bond, mag- 
netic and/or orbital order, and which of them is domi- 
nant depends on the details of the model relevant to the 
specific class of materials. 

In the present work we study a new mechanism for 
magnetic order in the S=l/2 pyrochlore lattice, driven by 
the Dzyaloshinsky-Moriya (DM) interactions. In the py- 
rochlore such interactions are expected to be present by 
symmetry. For the S=l/2 Heisenberg model on the py- 
rochlore lattice it has been suggested^ that the ground 
state is dimerized (non-magnetic), but macroscopic de- 
generacy still remains. For certain other lattices, such 
as the 2D pyrochlore and related models^ the ground 
state is a unique valence bond solid, and while the DM 
interactions (if present) can lead to non-trivial order in 
the ground state, such DM induced order can only occur 
above a critical threshold, due to its inherent competi- 
tion with the underlying dimer order£ In this work we 
show that in the 3D pyrochlore antiferromagnet, where a 
macroscopic degeneracy is present, the DM interactions 
have a more profound effect and can lift the degeneracy, 
leading to a chiral antiferromagnetic state with a small 
staggered magnetic moment. In an external magnetic 



field quantum transitions between weakly ordered states 
with different symmetries, depending on the field direc- 
tion, are possible. We determine the field-induced pat- 
terns for several field orientations, generally pointing in 
highly-symmetric crystal directions. The phase diagram 
at finite temperature is also briefly discussed. 




(a) (b) 

FIG. 1: (a) Pyrochlore lattice, (b) Distribution of DM vectors 
on a single tetrahedron (four of the six shown, see text). 

The spin Hamiltonian (S=l/2) is 

H = V]j ij S,S j + ^D ij .(S i xS j ), (1) 

where Dy are the DM vectors, to be specified later. 
We start by summarizing the results for Dy = 0, i.e. 
the Heisenberg case. Our starting point is the strong- 
coupling approach, similar to that of Refs. 00, with the 
lattice divided into two interpenetrating sub-lattices, one 
of them formed by "strong" tetrahedra (with exchange 
J), connected by "weak" tetrahedra (exchange J'). The 
"strong" tetrahedra then form a fee lattice, as shown in 
Fig. EJa), where every site represents a tetrahedron, and 
one can attempt to analyze the structure of the ground 
state starting from the limit J' <C J- 

For J' = the tetrahedra are disconnected, and on 
a single tetrahedron the ground state is a singlet and is 
twofold degenerate. We choose the two ground states as: 
ki> = -U[l,2][3,4] + [2,3][4,1]}, 1*2) = {[1,2][3,4]- 
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[2, 3] [4,1]}, where [k,l] denotes a singlet formed by the 
nearest- neighbor spins k and I, labeled as in Fig.^b). In 
the pseudo-spin T = 1/2 representation, so that T z = 1/2 
corresponds to \s\) and T z — —1/2 corresponds to \s2), 
one finds that third order is the lowest one contributing to 
the effective inter-tetrahedron Hamiltonian in the singlet 
sub-spaceA 




In the two-body part we have defined 

^03 = ^12 = 1/2, ^03 = ^12 = 1/ 6, 

^23 = ^01 = ^02 — ^13 = 1/3; 

n% = ngf = = -ng = 1/(2 Vz). (5) 

All remaining fiy = 0, i < j. The site indexes 
i,j refer to the fee lattice made of individual tetra- 
hedra, Fig. |^a), and it is sufficient to know the in- 
teractions on one " supertetrahedron" , shown in green 
(containing the sites 0,1,2,3). In the three-body in- 
teraction the indexes run over the values: (i, j , k) = 
{(3,2,1), (1,0,3), (2, 3,0), (0,1, 2)}. 

On a mean- field level the ground state of Ti- e f[ is defined 
by the following averages: 

(Tf ) = -V3/4, (7T) = 1/4, 

(Tf) = V3/4, (I?) = 1/4, 

(T$) = 0, (T|) = -1/2, 

(Iff) = 0, (T«f ) = 0. (6) 

This means that while a dimerization pattern sets in on 
sites 1,2,3, the pseudospins on the "0" sites, shown in blue 
in Fig. ffi a) remain "free" , i.e. there is no fixed dimer 
pattern on those sites and consequently a macroscopic 
degeneracy remains^ 

One should certainly keep in mind that the strong- 
coupling approach breaks artificially the lattice symme- 
try and while one hopes that the structure of the ground 
state is correct even in the isotropic limit J' = J, it is 
very difficult to assess this by other means (e.g. exact 
diagonalizations) at the present time. Nevertheless this 
approach is expected to provide reliable description of the 
ground state properties as long as the relevant physics re- 
mains in the singlet subspace, i.e. the triplet modes stay 
high in energy and no magnetic order is generated, as 
might be the case for the pyrochlore antiferromagnet due 



to the strong frustration. Fluctuations around the mean- 
field solution, Eq. (6), can lift the degeneracy, leading to 
unique dimer order. However the corresponding degen- 
eracy lifting energy scale is very small^ of the order of 
10 _3 /3, ee J' 3 /(48J 2 ). A unique (singlet) ground state 
is also produced if one starts the expansion from a larger 
cluster of 16 sites, with an ordering energy scale (energy 
gain) of 10 _2 J, extrapolated to the limit where all cou- 
plings are equals 

In what follows we will take the mean-field solution 
as a starting point and discuss a physical mechanism, 
based on the presence of interactions beyond Heisenberg 
exchange, that can lead to the lifting of degeneracy and 
consequently to (magnetic) order in the ground state. 




(a) (b) 

FIG. 2: (Color online.) (a) Fee lattice of tetrahedra (tetrahe- 
dron = dot) with interactions J' between them, (b) Antifcr- 
romagnetic chiral order on the blue (dark gray) tetrahedra, 
with magnetic moment |(S}| ~ D, induced by the DM inter- 
actions. On the gray tetrahedra (labeled as 1,2,3) the order 
has the same symmetry, but is much weaker |(S)| ~ D 3 <C D, 
Eq. (9), and is not shown. 



II. CHIRAL ANTIFERROMAGNETIC ORDER 
INDUCED BY THE DZYALOSHINSKY-MORIYA 
INTERACTIONS 

Now we consider the effect of the DM interactionsii*^ 
on the ground state properties. On a single tetra- 
hedron the DM vectors are distributed as shown in 
Fig. or explicitly: D 13 = ^(-1,1,0), D 24 = 

^(-1,-1,0), D 43 = ^(0,-1,1), D 12 = J|(0,-1,-1), 

D 14 = ^(1,0,1), D 23 = ^(1,0,-1). Here D is the 

magnitude of the (all equal) DM vectors. The directions 
of the DM vectors respect the pyrochlore lattice sym- 
metry and thus the DM interactions are expected to be 
always present in the system^ Since originate from 
the spin-orbit coupling; 11 ! 12 we have D -C J, J', and typi- 
cally the values of the DM interactions are several percent 
of the Heisenberg couplings. There are two DM distri- 
bution patterns that are equally acceptable on symme- 
try grounds - the one shown in Fig. ^b), and one with 
all directions of the vectors reversed (D^- — > -Dy). 
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These two cases were named, respectively, " indirect" and 
" direct" in Ref . The reader is referred to that paper 
for more details on Moriya's rules as applied to the py- 
rochlore lattice. In the extreme quantum case of S=l/2, 
and within our approach, we have found that the two 
allowed (by symmetry) DM distributions lead to qualita- 
tively the same physics (see discussion following Eq. (8)). 

Following the strong-coupling approach outlined above 
for the purely Heisenberg case, we have to determine how 
the singlet ground states |si), |s 2 ) on a single tetrahedron 
are modified by the presence of D. Since the DM in- 
teractions break the spin rotational invariance, they ad- 
mix triplets to the two ground states, not lifting their 
degeneracy^ We will also be interested in effects in the 
presence of an external magnetic field, and in this case 
the field (in combination with the DM interactions) also 
mixes certain triplet states with |si),|s2)- In order to 
determine the additional contributions to Ti- e Q, it is con- 
venient to express the spin operators on a singlet tetrahe- 
dron, labeled as in Fig. ^b) , in terms of the pseudospin 
operators. For magnetic field H = ^(1, 1,0) (along the 
1—3 bond), assuming D < J and H <C J, we obtain 
(defining the rescaled quantities D, H along the way): 

H = ^=(M,0); D = D/J, H = H/J, 
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The notation S,y simply combines in one line the for- 
mulas for both Si and Sj, where i and j label sites on 
a tetrahedron (as defined in Fig.^b)), while the left in- 
dex in Si j corresponds to the upper sign on the right 
hand side, and the right index - to the lower sign. The 
formulas Eq. (7) are obtained by using the ground state 
wave-functions, written explicitly in Ref. |8j (Eqs. (2,5)), 
to lowest order in D and DH. For magnetic field in the 
z direction, the corresponding expressions are given in 
Appendix A. 

First we analyze the case of zero magnetic field (H = 
0). Taking into account the connections between the 
tetrahedra (green bonds in Fig. EJa)), and using Eq. Q), 
we obtain an additional interaction term, so that the full 



effective Hamiltonian Tv 1 ^' becomes 



eff 



^ff M) =^cff-^^E^ 

i<j 



(8) 



where H e Q is the part originating from the Heisenberg 
exchanges, Eq. (2). The above result is obtained in low- 
est, first order in J' . While extra terms of the same power 
J'(D' j J) 2 also arise from the DM interactions D' on the 



inter-tetrahedral bonds, we find that they only give a 
small renormalization of the energy scale J' 3 /(8J 2 ) in 
Eq. (2) and are, therefore, neglected. 

We would like to also point out that in general the 
coupling constant in 7t! is not symmetric under 
D — > — D. To verify this requires a calculation of the next 
to leading order in D in Eq. (7). We have found that, in 
the case of zero field H = 0, the next order present is D 2 , 
and one has to substitute in all formulas D — > D %=D 2 . 

Consequently the same substitution has to be made in the 
coefficient — J'Z) 2 | in Eq. ©. Most importantly however 
the T^TV structure of the interaction is not affected by 

increasing the strength of D -C 1, and from now on we 
will work with the leading order in D. Therefore the 
physics (ground state structure) associated with the two 
DM distribution patterns will be the same. This conclu- 
sion might be connected with the fact that we have kept 
only the lowest non-trivial order in the coupling J' in the 
effective Hamiltonian - we have used this as our guiding 
principle as the difficulties associated with the derivation 
and analysis of higher orders seem insurmountable. 

We have performed mean-field calculations of the 
Hamiltonian defined by Eqs. f2 l3l4l8[) in the unit cell of 
Fig. |2{a), as represented by the four sites connected by 
green lines. The results can be particularly simply sum- 
marized in the limit D -C 1, which is also the case of 
physical relevance. It is physically clear that ferromag- 
netic order in the T? component is generated on the "0" 
sites, since no order in the T? ,z components (dimer or- 
der) was present on those sites without DM interactions 
(on mean-field level), Eq. (6). Indeed we find (Tq ) = 1/2, 
while for the other sites we have, to lowest non-trivial or- 
der in D, (Tf) « 1.8(L>/J') 2 ,i = 1,2,3. From Eq. © it 
is then clear that a non-zero average of the operator Tl y 
corresponds to a finite moment in the ground state, with 
mag nitude |(S) ; | = iV^lf). To summarize: 



PS) 

l<S)i| 



1/2, (2?>«1.8^j,i 



D 
71 



,i = 0; |(S)i| 



3.6 



1,2,3 
D 2 



^D — ,i = l,2,3. (9) 



Here |(S)i| stands for the magnitude of the moment on 
each site of pyrochlore lattice, belonging to a tetrahe- 
dron labeled by the index i. From (0) it follows that 
the moments point out of the cube's center (the cube is 
defined in Fig.^b)), leading to formation of sublattices 
and the order shown in Fig. |2fb) ■ Since from Eq. © 
|(S)i|/|(S) | - (D/J') 2 « 1, t = 1,2,3, we have ne- 
glected the magnetic order on those tetrahedra. 

The antiferromagnetic order of Fig. |2Ib) corresponds 
to non-zero scalar chirality \ — (S m • (S„ x S|)) 7^ 0, 
where to, n, I are any three spins on a given tetrahedron. 
The Ising symmetry T ; y — > is spontaneously bro- 

ken in the ground state, which in terms of real spins 
corresponds to the time-reversal symmetry broken state 
of Fig. |2b). In this state the two ground state wave 
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functions |<E>) and |^) (see (Al)) form linear combina- 
tions in the "crural" sector: a|$) + i/3\*if), where a, f3 
are real coefficients (a 2 + 2 = 1). This combination 
is ferromagnetically repeated on every T v ordered tetra- 
hedron. A straightforward calculation shows that both 
(Tj 7 ) ex a/3, x a P- The energy gain (per site of 
Fig. |2Ja)) from the formation of the ordered state is 
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(H { e g M) ) - (H eS ) « -1.8J'D 2 (D/J') 2 . The order 
we have just discussed is in competition with other mech- 
anisms for lifting of the degeneracy that could originate 
from the Heisenberg interactions themselves (e.g. fluctu- 
ations beyond the mean-field), typically also leading to 
very small energy scales^ 



III. MAGNETIC ORDER INDUCED BY 
EXTERNAL MAGNETIC FIELDS IN THE 
PRESENCE OF DM INTERACTIONS 

In the presence of an external magnetic field other pos- 
sibilities for lifting of the degeneracy exist. We will con- 
sider three symmetric field orientations, for which the 
results are particularly transparent. The magnetic field 
generally leads to splitting of the ground states, which in 
the language of the pseudospins produces an on-site "ef- 
fective magnetic field" h in the pseudospin z direction. 
The effective Hamiltonian has the form 



H 



(H) 

off 



(H) 

eff 



(10) 



We consider fields in the (1, 1, 0) and (0, 0, 1) directions, 
as well as comment on the case (1,1,1), where the axes 
are defined in Fig. ^b). Using the wave- functions in a 
field we obtain (see Eq. (A3)): 



1D 2 H 2 /J*, H = $(1,1,0) 
-D 2 H 2 /J 3 , H = #(0,0,1) 
0, H= $(1,1,1) 



(11) 



in l|10|) represents lattice contributions, originat- 
ing from the various combinations in Eq. Q once the 
tetrahedra are coupled, and also producing terms of or- 
der D 2 H 2 . These terms are cumbersome and are not 
explicitly written, but their effect is taken into account 
in the (numerical) mean-field implementation within the 
unit cell of Fig. Ufa). A further discussion appears in 
Appendix B. 

We will mostly discuss the two cases with h ^ 0. Then 
the on-site h term in Eq. (|10fl is responsible for the main 
effect, namely competition between order in the T? pseu- 
dospin component and order in the "crural" T? compo- 
nent favored by Eq. JHJ. Therefore the physics is that 
of the transverse field Ising model (although in our case 
the unit cell is larger) . It is also clear that the mentioned 
competition is most effective on the "0" (blue) sites, while 
the non-zero averages of T. ' on the other sites are not 
much affected by the presence of small D and H. We have 





H=0 



0<H<H„ 



H>H„ 




FIG. 3: (Color online.) (a,b,c) Evolution of magnetic order 
on the blue (dark gray) tetrahedra of Fig. 2(b) in an external 
magnetic field in the (1,1,0) direction. Upper row: field- 
induced order on the rest of the tetrahedra. The tetrahedra 
are labeled 0,1,2,3 as in Fig. 2(a,b). Blue (dark gray) arrows 
in (a,b) correspond to moments | (S) | ~ D, while the red (light 
gray) arrows on the upper row and (b,c) are the field-induced 
moments |(S)| ~ DH. 



found that a quantum transition takes place between a 
state with (T<f) ^ 0,H < H c and (T a ) = 0,H > H c . 
The result for D <C 1 can be written in an explicit way, 
and we have for the field H = $(1, 1, 0) 



<7dV = 7 




(T o r = l/4-(T ^) 



H <H r 



y\2 



5.3 A /-£> (12) 



(13) 



(and (Tq) < since h > 0). The values of the spin 
moments for given values of (T x,v,z ) on a tetrahedron 
can be determined directly from Eq. J7J. On the "0" 
(blue) sites this leads to evolution of the magnetic order 
as shown in Fig. |3Ja,b,c). For H = there is only chiral 
order (blue arrows) with moment |(S)| ~ D, changing, 
for H > into a combination of chiral and field induced 
order (red arrows) with | (S) | ~ DH. Gradually, as H ap- 
proaches H c the chiral order diminishes (Eq. leav- 
ing for H > H c only the field- induced component, equal 
to | <S> | = Ml(T *)| = DH/2, H > H c . 

On the tetrahedra 1, 2, 3 labeled as in Fig.[2Ia,b) there 
is virtually no evolution as a function of the field, and 
the order is determined by Eq. J7J with (Ti X ' z ) fixed by 
the Heisenberg exchanges, see Eq. (6). This leads to the 
magnetic moments (proportional to DH) shown in Fig. [3 
upper row. On tetrahedra 1 and 2 the spins point along 
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the internal diagonals of the cube perpendicular to the 
field. Dimerization is also present in the ground state 
(bolder lines = stronger bonds) and co-exists with the 
magnetic order. 



H 




H=0 0<H<H p ,, H>H C 

C U u 



10,0,1J 




FIG. 4: (Color online.) Same as Fig. 3, but for an ex- 
ternal magnetic field in the (0, 0, 1) direction. Blue (dark 
gray) arrows in (a) correspond to the DM induced order with 
|(S)| ~ D, red (light gray) arrows on the upper row and (c) 
correspond to the field- induced component |(S)| ~ DH, and 
black arrows in (b) are a mixture of the two. 

A similar quantum transition takes place for magnetic 
field in the z direction H = H(0, 0, 1). In this case the 
formulas (A2) from Appendix A have to be used, and the 
field-induced order is shown in Fig.0] The critical field is 
also somewhat smaller in this case H c ~ 3.8( J/ J') 1/,2 D, 
mainly due to the fact that h is larger by a factor of 2, see 
Eq. (11). For other, less symmetric field directions, the 
form of the effective Hamiltonian, and consequently the 
field-induced patterns can be quite complex. Finally, in 
the case of a field H ~ (1, 1, 1), when h — 0, the quantum 
transition described above does not take place, and the 
chiral order of Fig. OJa) essentially does not evolve. In 
this case the various additional terms similar to the ones 
described in Appendix B may lead to small, sub-leading 
deviations from the perfect chiral sate. 

In addition to the field-induced ordered patterns of 
Fig. |21 and Fig. 0] determined mostly by the inter- 
tetrahedral interactions, a single tetrahedron with DM 
interactions also possesses a finite moment in the direc- 
tion of the field^ meaning that the spins in Fig. [3] and 
Fig. 2| would also tend to tilt in that direction. However 
the moment along the field is proportional to D 2 H, as can 
be deduced from the fact that the ground state energy 
varies as D 2 H 2 from (A3). Consequently this component 
has not been taken into account in Eqs. 0A2), valid to 
lowest order in D, H. Finally, we emphasize that while 
we have assumed D,H to be small, the ratio D/H can 



be arbitrary, meaning that the quantum transitions in a 
field are within the limit of validity of our approach. 



IV. PHASE DIAGRAM AND DISCUSSION 

At finite temperature we expect the phase diagram to 
look as presented in Fig. [5] (it is assumed that h ^ 0). The 
higher transition temperature T c \ ~ J' 3 / J 2 corresponds 
to the scale below which the translational symmetry is 
broken (dimerization occurs), and is determined by the 
energy scale in Ti. e Q for D = 0, Eq. (2). We expect T cl 
to have weak dependence on magnetic field when DM 
interactions are present. At a lower scale T c2 the Ising 
T v — > — T y symmetry is spontaneously broken by Eq. JSJ. 
For H = we can estimate T c2 ~ J'D 2 (D/J') 2 . At 
fixed field this finite-temperature transition is in the 3D 
Ising universality class, and the specific heat diverges as 
C ~ |T - T c2 (H)\- a , a re O.llM We emphasize that 
Fig. [5] shows only the low-field part of the phase diagram 
(since H c ~ D -C J), while the physics at high fields 
cannot be determined within the effective Hamiltonian 
framework presented here. 

ii T 

<T X ' Z > = 



T 

_c_1 

<T X ' Z > 4 partial dimer order 




FIG. 5: Schematic phase diagram at non-zero temperature 
in the presence of small magnetic field (H <C J) and DM 
interactions. 

In certain pyrochlores, such as the gadolinium tita- 
nium oxides with S=7/2, field-driven phase transitions 
have been observed^ although in this material the mag- 
netic order is typically explained as originating from the 
long-range dipolar interactions. For such large value of 
the spin the DM mechanism for magnetic order, at least 
the way it is developed in this work, should not be effec- 
tive since our calculations were based on strong singlet 
correlations in the ground state. At the moment it is hard 
to point out a class of materials where the DM interac- 
tions are definitely expected to be dominant with respect 
to other anisotropics capable to produce ordering; some 
possible examples are given in Ref. Il3l In particular, 
our results are specific to the case S = 1/2, while most 
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currently known pyrochlores have higher S. Let us also 
point out the main differences between the present work 
and the purely classical models (a.) We have found 
that the antiferromagnetic order is weak, determined by 
the DM interaction scale itself, (b.) The "chiral" non- 
coplanar pattern is the stable one for S = 1/2 in the 
absence of a field. Field-induced patterns then dominate 
for H > H c ~ D. Let us also give some estimates: if we 
take the optimistic viewpoint and apply our formulas to 
the case J = J', and take J ~ 100K,D/J ~ 10 _1 , then 
the characteristic temperature for onset of chiral order 
would be T C 2 ~ lOmK, and the magnitude of the mo- 
ment |(S)| ~ 0.1. The characteristic critical fields would 
be H c ~ lO-fsT. These should be viewed as order of magni- 
tude estimates. Inevitably the scale T C 2 falls into the mK 
range, which, in combination with the smallness of (S) it- 
self, would probably make the chiral state rather hard to 
observe with neutron scattering techniques used to probe 
the spin structured The spectral weight of magnetic ex- 
citations in such neutron scattering measurements would 
be determined by D itself. However the scale of H c of 
around 10T or smaller (also determined by D), means 
that the various field-induced patterns, strongly depen- 
dent on the magnetic field direction, might be accessible. 

In conclusion, we have shown that DM interactions 
can induce weak antiferromagnetic order characterized 
by non-zero chirality. In an external magnetic field quan- 
tum transitions between the chiral state and field-induced 
ordered states take place. Field-induced patterns with 
different symmetries depending on the direction of the 
field are very characteristic of the presence of DM inter- 
actions. We have used an expansion around a configu- 
ration which breaks the lattice symmetry^ and leaves a 
macroscopic degeneracy, subsequently lifted by the DM 
interactions. Full restoration of lattice symmetry within 
such an approach seems impossible to achieve, as it is im- 
possible for example in the large-N approach^ Neverthe- 
less we expect that without DM interactions the ground 
state properties and the inherent degeneracy present in 
this strongly frustrated magnet are well accounted for. In 
this situation the DM interactions "push" the pyrochlore 
lattice towards the ordered states analyzed in the present 
work. More generally, the DM interactions can be rele- 
vant and lead to weak antiferromagnetism in strongly 
frustrated systems, where the Heisenberg exchanges on 
their own fail to produce long-range order. We also em- 
phasize that the physics behind the weakly antiferromag- 
netic states with different symmetries discussed in this 
work is very different from the phenomenon of weak fer- 
romagnetism, usually associated with DM interactions. 



APPENDIX A: SPIN OPERATORS AND 
ENERGY LEVELS IN MAGNETIC FIELD 

Here we present the expressions for the spin operators 
for magnetic field in the z direction. The two ground 
state wave-functions |si), \s\) are modified in the follow- 
ing way in the presence of DM interactions (\si) — > |$), 

l*a>-l*»: 

H = #(0,0,1); 

!•>- W + ^IW-W + W + I*>] 
Vqdh 

+ ^ UPx) + \Py) - \lx) + \q y )\ , 

l*> = \s2) + ^[\p x ) + \p v ) + \q x )-\q v )]+W\t z ) 

J J J! 

+ ^ HP*> + \Pv) + Ife) + \ly)) > ( A1 ) 

where jp^), \q/x), \t^), /i = x,y,z are the three excited 
triplet states on a tetrahedron. From these equations we 
obtain 
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QX 

D 2,4 






DH 

V6 


QV 
°2,4 


\/6 




DH 


QZ 

J 2,4 


- _ 2 A T y 
V6 ' 







(A2) 



We also give the ground state energy splitting on a 
single tetrahedron, in (weak) external magnetic field with 
arbitrary direction 



H — (H x , H y , H z ) 

3J D 2 |H| 2 
~2 2J~ 



E, 



(1,2) 



J 3 



D' 2 



± 2J^ V |H|4 " * {H * H v + H * H * + H " H y ] 

(A3) 

From here the "effective magnetic field" h appearing in 
the pseudospin Hamiltonian (10) is: h = E^ — Eq 1 ^. 
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APPENDIX B: EFFECTIVE HAMILTONIAN 
CONTRIBUTIONS IN MAGNETIC FIELD 



We briefly discuss the structure and treatment of the 
term a in (10). As is clear from (7) and (A2), po- 



7 



tentially contributions of two types appear: (a.) terms 
of order D 2 H, and (b.) terms of order D 2 H 2 . While 
treating these terms we will assume that the unit cell 
structure of Fig. 2(a) does not change. Indeed, since we 
are interested in the case of weak fields H <C 1, the above 
coupling constants are of sub-leading order with respect 
to the case without a field, and therefore one does not 
expect the unit cell to change. Under this assumption 
we find that the terms of order D 2 H vanish identically 
for the field directions considered in this work. The re- 
maining contribution, e.g. for a field H = -^=(1,1,0), 

written (per site) within the unit cell of Fig. 2(a) is 



where 



6 



D 2 H 2 6H, 



(Bl) 



SH = -4(T$ + Tg)(T? + Tg)- 
v/3(T *-Tf)(Tf-T|)- 
y/3(Tg-Ti)(T?-T?)- 
3(Tgl% - T Z T 3 Z + TfTf 



TfT|). (B2) 



These terms were taken in to account in our numerical so- 
lution of the mean-field equations corresponding to (10). 
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